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ABSTRACT 

Aims, at Centauri (NGC 5139) contains large numbers of variable stars of different types and, in particular, more than a hundred 
RR Lyrae stars. A homogeneous (in terms of instrument, image quality, and time coverage) sample of high-quality near-infrared 
(NIR) RR Lyrae light curves could thus be gathered by performing an extensive time-series campaign aimed at this object. We have 
conducted a variability survey of cd Cen in the NIR, using ESO’s 4.1m Visible and Infrared Survey Telescope for Astronomy (VISTA). 
This is the first paper of a series describing our results. 

Methods. cd Cen was observed using VIRCAM mounted on VISTA. A total of 42 and 100 epochs in J and K s , respectively, were 
obtained, distributed over a total timespan of 352 days. Point-spread function photometry was performed using DAOPHOT and 
DoPhot in the inner and outer regions of the cluster, respectively. Periods of the known variable stars were improved when necessary 
using an ANOVA analysis. 

Results. An unprecedented homogeneous and complete NIR catalogue of RR Lyrae stars in the field of cd Cen was collected, allowing 
us to study, for the first time, all the RR Lyrae stars associated to the cluster, except 4 located far away from the cluster center. 
Membership status, subclassifications between RRab and RRc subtypes, periods, amplitudes, and mean magnitudes were derived for 
all the stars in our sample. Additionally, 4 new RR Lyrae stars were discovered, 2 of them with high probability of being cluster 
members. The distribution of a> Cen stars in the Bailey (period-amplitude) diagram is also discussed. Reference lines in this plane, 
for both Oosterhoff type I (Ool) and Oosterhoff type II (OoII) components, are provided, both in J and K$. 

Conclusions. In the present paper, we clarify the status of many (candidate) RR Lyrae stars that had been unclear in previous studies. 
This includes stars with anomalous positions in the color-magnitude diagram, uncertain periods or/and variability types, and possible 
field interlopers. We conclude that to Cen hosts a total of 88 RRab and 101 RRc stars, for a grand total of 189 likely members. We 
confirm that most RRab stars in the cluster appear to belong to an OoII component, as previously found using visual data. 

Key words. Stars: variables: RR Lyrae - Stars: variables: general - Stars: horizontal-branch - Galaxy: globular clusters: individual: 
a> Centauri (NGC 5139) - Infrared: stars - Surveys - Catalogs 


1. Introduction 


RR Lyrae stars (RRLs) are well known for being reliable dis¬ 
tance indicators and tracers of the old stellar populations of 
galaxies ( |Smith|1995| |Catelan|2009| |Bono et al.|201 1| and ref¬ 
erences therein). In the dawning era of time-series infrared (IR) 
photometric surveys, new observations in the near- and mid-IR 
have recently been providing large and homogeneous samples of 
RRLs in different environments (e.g., |Dekany et al.|2013||Klein| 
et al. 2014), which can be used to improve our understanding 
of this variability class and as tracers of structure in the nearby 
Universe. 


* Based on observations collected at the European Organisation for 
Astronomical Research in the Southern Hemisphere, Chile, with the 
VISTA telescope (project ID 087.D-0472, PI R. Angeloni). 


Especially noteworthy, in this context, is the fact that RRLs 
follow well-defined period-luminosity (PL) relations in the IR 
(e.g., Longmore et al. |1990| Catelan et alT[ |2004| . Furthermore, 
the interstellar extinction in Ks is approximately a tenth of the 
one in the optical, thus making it possible to observe RRLs even 
in the most obscured regions of our Galaxy, e.g., its inner bulge. 
Despite these advantages, there are still few extensive variability 
studies in the near-infrared (NIR), due to a number of factors, in¬ 
cluding the lack, until recently, of large-format detectors, which 
rendered NIR surveys extremely costly, in terms of observing 
time. In addition, the low amplitudes of many types of vari¬ 
ables, including RRLs, in the NIR, compared to the optical, have 
hindered the detection of new variables using these filters. The 
situation has only recently started to change, particularly with 
the start of operations of ESO’s 4.1m Visible and Infrared Sur- 
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vey Telescope for Astronomy (VISTA), located in Cerro Paranal, 
Chile ( |Emerson et al.|2006| [Emerson & Sutherland|2010| l. 

The VISTA Variables in the Via Lactea (VVV) ESO Public 
Survey is an extensive variability study in the AC,-band of the 


Galactic bulge and an adjacent portion of the disk (see Minniti 
et al.|2010[|Catelan et al.|2011||2013| and references therein). In 
order to automatically classify the large samples of variable stars 
detected by this survey, our team started to build a database of 
well-defined, high-quality NIR light curves for various variabil 
ity classes - the VVV Templates Project (Catelan et al.||2013 


Angeloni et al. |2014j >. In this context, to Centauri (NGC 5139) 


a comprehensive assessment of their number, types, membership 
probabilities, and light curves, and also improved light curves 
for several RRL variables. We point out some erroneous cross¬ 
identifications between the We I drake et al.| (|20Q7[> and |Kaluzny| 


et al. (2004) catalogues, and propose a change in classification 


was identified early on as one of the most suitable targets to ob¬ 
tain template light curves, given that it hosts a large number and 
variety of variable stars. 

to Cen stands out as one of the three most RRL-rich glob¬ 
ular clusters known, only surpassed by M3 (NGC 5272, with 
Clementini et aT||j004| and M62 (NGC 6626, with 
Contreras et al.|2005 2010). From the several vari- 


230 RRLs; 

224 RRLs;_ 

ability studies that have been carried out in the optical over the 
years with different telescopes and temporal baselines, a large 
(although still incomplete) catalogue of RRLs in to Cen has been 
produced (Clement et al. 2001 [ C. M. Clement’s online cata- 
iogue02Ol 3 edition, hereinafter C13). A wide-held, NIR vari¬ 
ability survey of to Cen allows us to derive and calibrate the 
NIR PL relations of several different types of stars, including 
SX Phoenicis, type II Cepheids, and RRL stars. In this sense, 
our study of to Cen with VIRCAM@ VISTA has provided the 
largest homogeneous NIR sample of RRLs ever collected in any 
globular cluster, leading to tight PL relations, as described in our 
next papers in this series on to Cen variables. 

a> Cen is the most massive and luminous globular cluster 
in the Milky Way, and it may have been even more massive in 
the past (e.g., |Valcarce & Catelan|2011( and references therein). 
Nowadays, to Cen is well known for hosting at least two stellar 
populations with metallicity peaks centered at [Fe/H] ~ -1.6 
and -1.1 dex ( |Bedin et al.||2004[ |Joo & Lee||2013| and ref¬ 
erences therein). In addition to presenting variations in [Fe/H] 
that may be correlated with age, it has been found that to Cen’s 
stars also present variations in the abundances of many other el¬ 
ements, both heavy and light, in addition to helium (at a level 
AT ~ 0.17 dex) and age (e.g., |Norris||2004HPiotto et al.||2005[ 
[Dupree et al.|20lT]|Valcarce & Catelan|2011|). 

Despite the fact that He enhancement appears to be required 
to explain the lower color-magnitude diagram (CMD) morphol¬ 
ogy of to Cen, an impact on the cluster’s variable stars has not yet 
been established. |Sollima et al .| ( f2006[ > found that the RRL stars 
with metal-intermediate composition ([Fe/H] ~ -1.2), fainter 
than the bulk of the dominant metal-poor population (with pri¬ 
mordial helium abundances and [Fe/H]-1.7 dex), are in good 

agreement with the corresponding horizontal branch (HB) mod¬ 
els with Y — 0.246. More recently, Marconi et al. (2011 1 studied 
the impact of the helium content on the RRL properties based on 
evolutionary and pulsation models, and again could not find evi¬ 
dence of He enhancement among the RRL stars. The most likely 
interpretation of these results is that any He-enhanced popula¬ 
tions in the cluster likely end up far to the blue on the zero-age 
HB, thus being unlikely to populate the classical instability strip 
(e.g.,[joo & Lee|2013[ ). 

In this paper, in addition to describing our VISTA observa¬ 
tions, we also update the census of RRLs in to Cen. In this sense, 
we provide, using our data and comparing with previous studies. 


for a few stars. The membership status of all the suspected field 
stars is revisited. Last but not least, we present four new RRLs, 
two of which we classify as genuine cluster members. In future 
papers of this series, we will study the NIR PL relations for RRL, 
type II Cepheids, and SX Phe stars, and present the full updated 
catalogue of variable stars (including also eclipsing binaries) in 
the to Cen field. 


2. Observations 

A total of 42 and 100 epochs in J and AT,, respectively, were 
taken using VIRCAM (|Dalton et al.||2006|) mo unted on VISTA 
( [Emerson et al.||2006[ [Emerson & Sutherland||2010l ). All of the 
images were centered on the cluster center. The heart of the 
VIRCAM camera is a 4 x 4 array of Raytheon VIRGO IR de¬ 
tectors (2048 x 2048 pixels), with a pixel size of 0734 (Dalton 
|et al.|2006] ). The size of a uniformly covered field (also called a 
“tile”) is 1.501 deg 1 2 , thus allowing unprecedented spatial cover¬ 
age in the NIR bandpasses. The effective field of view (FoV) of 
the camera is larger than the previous studies of Kaluzny et al. 
P004| ) and Weldra ke et al.| ( |2007) ), both carried out in the visi¬ 
ble, but not large enough to cover to Cen in its whole extension 
(its tidal radius is - 52'; e.g., Ferraro et al. 2006 and references 
therein). Our observations have a baseline of 352 days, which 
provides an unprecedented phase coverage in the NIR. 

The characteristics of the observations, data reduction and 


photometry extraction are the same as those explained in Navar- 
rete et al. ( |2013| >. In particular, the images were reduced by 


1 http://www.vvvtemplates.org 

2 http://www.astro.utoronto.ca/~cclement/read.html 


the Cambridge Astronomy Survey Unit (CASUjj^and crowded- 
field, point-spread function (PSF) photometry was performed us¬ 
ing the DAOPHOT II/ALLFRAME (|Stetson||19~87l [1994] ) and 
DoPhot ( |Schechter et al.|1993[ [Alonso-Garcla et al.|2012[ ) pho¬ 
tometry packages for the cluster central (i.e., the innermost ~ 
10') and outer regions, respectively. 

The DoPhot and ALLFRAME photometries were calibrated 
independently, using in both cases the aperture photometry pro¬ 
vided by CASU, which is in the VISTA photometric system. As 
described in detail in [Minniti et al.| ( |2010] >, the CASU pipeline 
determines the photometric zero points using 2MASS stars 
( jSkrutskie et al.|2006j ) in the pawprint frames, then transform¬ 
ing these 2MASS magnitudes into VISTA magnitudes according 
to the equations derived by Hodgkin et al. ( 2009[ ). Magnitudes 
for common stars derived using both packages were then com¬ 
pared, finding an excellent agreement in the magnitude calibra¬ 
tion (at the level of RRL stars, the average AT, -band magnitude 
difference between the DoPhot and ALLFRAME magnitudes is 
~ 0.01 mag only). 


3. Results 

Time-series photometry of most of the previously known vari¬ 
able stars listed by C13 (from VI to V410) and |Weldrake et ah 
( j2007| ) was recovered using a cross-match between those cata¬ 
logues and our complete PSF photometric catalogue (i.e., pho¬ 
tometry for all the stars at all the observed epochs). The match¬ 
ing radius was set at 1". When one star was not recovered under 
this procedure, the Samus et al.j ( 2009} catalogue was used as 


http://casu.ast.cam.ac.uk/surveys-proj ects/vista 
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reference. In the Samus et al.| catalogue, variable stars from VI 
to V293 are listed, most of them with updated coordinates based 
on 2MASS observations. The cross-match using the Samus et al. 


catalogues allowed us to recover 16 variable stars (including 8 
RRL, 3 SX Phoenicis, 4 eclipsing binaries, and 1 spotted vari¬ 
able) that were not recovered using the C13 catalogue. When 
neither of these reference catalogues gave a match for a cer¬ 
tain variable star, an increasing matching radius was used, un¬ 
til a known variable star was found. To pinpoint the latter, the 
time series for all stars that were inside the matching radius 
were phased using the star’s period as listed by Cl3, until rea¬ 
sonable agreement was found. As a consistency check, periods 
for those stars were also calculated using the ANOVA algorithm 


(Schwarzenberg-Czerny 1989), in an effort to determine if the 


known variable is present in the area covered by the matching 
radius but with a different period according to our observations. 
This approach led to the discovery of V411 as a new RRL star 
(Navarrete et a l.|2013[ ), which Weldrake et al.| ( 2007j l had previ- 
ously (but incorrectly) claimed to be the same as VI44 from the 
Kaluzny et al. ( |2004j > catalogue. 

Figure |T] shows the to Cen Ks, (J - Ks) CMD, as ob¬ 
tained from our ALLFRAME photometry for the central regions 
(see Sect. [2|. In this diagram, RRLs from the whole observed 
area are superimposed as red and blue circles for RRab (i.e., 
fundamental-mode) and RRc (first-overtone) types, respectively. 
Star symbols indicate variables from |Weldrake et al.| ( |2007| 
(see Sect. 0 - and new variables are shown as open circles (see 
Sect. |3.3| l. Most of the RRLs are located in the horizontal branch 
(HB), as expected for bona-fide cluster members, but a few stars 
fall close to the main sequence, or even at the turn-off point. A 
zoom-in around the RRL instability strip (IS) region is shown in 
the bottom panel of Figure [I] where there are also two type II 
Cepheids, marked as green triangles. All the RRL stars labeled 
in both figures were considered in our analysis and studied in 
terms of their periods, variability types, and cluster membership. 
In the following, we clarify the status of many of the RRL stars in 
to Cen that had been unclear in previous studies. This includes 
stars with anomalous CMD positions, uncertain periods or/and 
variability types, and possible field interlopers. 

Table [TJlists the characteristics of 21 previously known RRL 
that either have anomalous positions in the CMD, were sug¬ 
gested as field stars by previous studies (and fall within our 
FoV), or do not have periods (or even variability type) listed in 
the literature. Column 1 lists the ID of the variables according 
to C13, except for the two stars from Weldrake et al.|(2007} ), be¬ 
cause they are not included in the C13 catalogue (ID-91, ID-99); 
columns 2 and 3 show the coordinates based on VISTA astrom¬ 
etry (which in turn is tied to the 2MASS system); columns 4 and 
5 display the intensity-weighted mean J and Ks magnitudes, re¬ 
spectively; column 6 gives our derived periods, while the vari¬ 
ability type (RRab, RRc, eclipsing binary [EB], or anomalous 
Cepheid [ACEP]) is given in column 7; the distance to the clus¬ 
ter center for each variable is listed in column 8. The member¬ 
ship probability for the stars, according to the two main proper 
motion studies of to Cen, are listed in columns 9 and 10, as pro¬ 
vided by Ivan Leeuwen et al. (2000 hereafter vLOO) and |Bellim| 


et al. 


(2000 hereafter B09), respectively]^] Finally, column 11 


van Leeuwen et al. 


([2000) constructed a proper motion catalogue 


for stars in to Cen based on observations of photographic plates of the 
cluster, taken between 1931 and 1983 with the Yale-Columbia 66-cm 
telescope. The covered area has a radius of ~ 30' (about 60% of the 
cluster’s tidal radius). The |Bellini et al.|(2009fr proper motion catalogue 
was in turn derived from astrometry based on CCD images taken with 
the Wide-Field Imager on the 2.2-m telescope at ESO La Silla, with a 


14 


. RRc 

• RRab 
o NV 

* W07 


V143 

v5 - 2 * *, NV455 



- 0.2 0.0 0.2 0.4 0.6 0.8 1.0 

J-K, 



J-K, 

Fig. 1. Top: CMD of the central part of the cluster (11 arcmin 2 , grey 
dots), ab- and c-type RRL are marked as red and blue circles, respec¬ 
tively. RRL from Weldrake et al. (2007) are marked as star symbols. 
Open circles mark the new RRLs found in this work, with NV455 and 
NV456 (labeled) likely being cluster members (see Sect. |3.3| i. Field, 
blended, or erroneously classified stars are labeled and explained in the 
text. Bottom: Zoom-in around the RRL region in the CMD of to Cen, 
with a few noteworthy stars labeled. V43 and V60, marked with green 
triangles, are known type II Cepheids. 


shows the membership status according to our photometric study 
(m: member, f: field). A brief description of each of these vari¬ 
ables is presented in the Appendix. 

precision of ~ 7 mas. The B09 observations were carried out with a 
baseline of 4 years, and extend to a limiting magnitude of B ~ 20 mag, 
which is up to four magnitudes deeper than previous studies. 
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Table 1 . Properties of member and field RRLs with previously uncertain status 


ID 

RA 

(J2000) 

DEC 

(J2000) 

(J) 

(mag) 

(K s ) 

(mag) 

p(a) 

(days) 

Type 

d< b > 

(arcmin) 

rpv LOO 

r H 
(%) 

®B09 

A* 

(%) 

Memb. 

status 

V19 

13:27:30.12 

-47:28:05.74 

13.852 

13.631 

0.299551 

RRc 

7.27 

100 

100 

m 

V21 

13:26:11.15 

-47:25:59.30 

13.538 

13.362 

0.380812 

RRc 

6.71 

100 

99 

m 

V52 

13:26:35.15 

-47:28:04.33 

12.831 

12.634 

0.660386 

RRab 

2.16 

53 

45 

m 

V56 

13:25:55.44 

-47:37:44.44 

13.612 

13.301 

0.568023 

RRab 

12.53 

98 

100 

m 

V68 

13:26:12.79 

-47:19:36.12 

13.168 

12.897 

0.534696 

RRc/ACEP? 

10.87 

100 

100 

m? 

vso"^ 

13:28:54.96 

-47:30:16.42 

13.568 

13.360 

0.37718 

RRc 

21.62 

out 

out 

m 

V84 

13:24:47.47 

-47:49:56.18 

13.031 

12.712 

0.579873 

RRab/ACEP? 

29.24 

0 

out 

m? 

V143 

13:26:42.59 

-47:27:28.98 

12.783 

12.510 

0.820734 

RRab 

1.51 

90 

99 

m 

V151 

13:28:25.31 

-47:16:00.01 

13.455 

13.211 

0.407756 

RRc 

20.94 

out 

out 

m 

V168 

13:25:52.74 

-47:32:03.19 

14.164 

13.916 

0.321299 

RRc 

9.78 

0 

0 

f 

V172 (d> 

13:27:55.04 

-47:04:38.50 

13.220 

12.927 

0.7380 

RRab 

26.72 

out 

out 

m 

V173 

13:29:43.13 

-47:16:53.80 

13.624 

13.360 

0.3590 

RRc 

32.05 

out 

out 

m 

V177 w 

13:29:04.15 

-47:36:21.42 

13.690 

13.479 

0.31469 

RRc 

24.31 

out 

out 

m 

V180 

13:28:15.22 

-47:40:32.02 

13.201 

12.839 

0.39006 

WUMa 

18.93 

0 

out 

f 

V181 

13:30:00.35 

-47:48:44.96 

14.823 

14.535 

0.5884 

RRab 

38.17 

out 

out 

f 

V183 

13:29:39.50 

-47:30:18.36 

14.690 

14.449 

0.29603 

RRc 

29.13 

out 

out 

f 

V268 

13:26:35.11 

-47:26:11.15 

13.156 

12.827 

0.812922 

RRab 

3.30 

70 

99 

m 

V283 

13:27:36.53 

-47:46:40.44 

16.947 

16.682 

0.517233 

RRab 

19.73 

out 

out 

f 

V411W 

13:26:40.77 

-47:28:17.00 

13.123 

12.796 

0.8449 

RRab 

1.20 

100 

99 

m 

V433^ 

13:29:03.59 

-47:48:58.25 

14.737 

14.151 

0.6681 

? 

30.58 

out 

out 

f 

ID-91 

13:26:14.03 

-47:53:08.20 

15.782 

15.430 

0.637238 

RRab? 

25.00 

out 

out 

f 

ID-99 

13:24:40.44 

-47:45:22.54 

15.777 

15.275 

0.670361 

? 

27.07 

out 

out 

f 


Notes. (fl) Periods are from Kaluzny et al. (20041, except for V80, V172, V173, V177, V180, V181, V183 and V411, which were derived by us. 
,h> Distance to the cluster center was calculated usin g a> Cerfs center as rep orted by SIMBAD: RA = 13:26:47.28, DEC = -47:28:46.1, both in 
J2000. (c) ID-53 in |Weldrake etall(|2007|>. w ID-18 injWeldrake et al.||2007| i. (e) ID-49 in |Weldrake etaL| ( |2007 | >. ,f) V41 1 has the same coordinates 
a nd periods as first repor ted by Navarrete et al.|j20 1 3j. Mean magnitudes differ from that paper becau se the values presented in this paper, but not 


Navarrete et al. 1 2013 1 , are intensrty-averaged mean magnitudes. (s) ID-74 in Weldrake et al. 1 2007 1 . 


3.1. Bailey diagram 


As shown by Clement & Rowe (2000) based on visual data, indi¬ 
vidual stars in u> Cen predominantly follow a locus on the Bailey 
diagram (period-amplitude plane) that is characteristic of Oost- 
erhoff type II (OoII) systems, with individual RRL typically hav¬ 
ing longer periods at a given amplitude than their Oosterhoff type 
I (Ool) counterparts. The details, systematics, astrophysical in¬ 
terpretation and importance of the Oosterhoff phenomenon are 
extensively discussed in Catelan|( [2009] l. 

Figure [2] shows, for the first time, the distribution of oj Cen 
RRL stars in the NIR Bailey diagram, both in the J (top panel) 
and in the Ks (bottom panel) bands. RRab stars with documented 
Blazhko effect (Blazhko 1907) in the catalogue of Kal uzny et al.| 
(2004) are plotted as open red circles. Several features can be 
apprehended from these plots. First, the positions occupied by c- 
and ab-type RRL stars are clearly very different, as also seen in 
the visual. Second, there is significantly more scatter in the K$- 
band Bailey diagram than is present in J. This is, at least in part, 
due to the smaller amplitudes in K$, compared to J. Third, and 
as it is also commonly seen in the visible, the Blazhko stars tend 
to have smaller amplitudes, at a given period, than their non- 
Blazhko counterparts. This is particularly noteworthy as the ac¬ 
tual amplitude modulation that is brought about by the Blazhko 
effect is not clearly seen in our NIR data for any of u> Cen’s RRab 
stars, which could be either because the amplitude changes are 
too small in the NIR (notwithstanding the fact that Fig.[2]might 
hint at amplitude modulations exceeding 0.1 mag in at least some 
cases), and/or because our time coverage is insufficient to re¬ 
veal long-term modulation (recall that our observations have a 
timespan of ~ 300 days, with gaps). Indeed, to the best of our 


knowledge, there have been no clear-cut measurements of the 
Blazhko effect in the IR, even though [Sollima et ak] ( j2008[ ) have 
reported increased scatter in their NIR light curves of RR Lyr it¬ 
self, which they attribute to the star’s well-known Blazhko-type 
light curve modulation. Fourth, there is a tendency for the locus 
occupied by ab-type RRL to flatten out, as one moves from bluer 
to redder wavelengths. In fact, this trend appears to continue to¬ 
wards the mid-IR regime, where the relation eventually becomes 
completely flat (Gavrilchenko et al. 2014). Fifth, and again as 
previously emphasized by |Clement & Rowe| ( {2000| i, most RRab 
stars in oj Cen tend to clump around a fairly well-defined locus, 
particularly in the /-band Bailey diagram. Based on our J am¬ 
plitudes for the “OoII-like” RRab stars, this locus is fitted by the 
following expression: 

= 0.064 - 2.481 log P+ 10.345 (logP) 3 , (1) 

where the period is given in days. The amplitude-period corre¬ 
lation given by Equation [I] is plotted as a dashed line in the /- 
band Bailey diagram (Figure [2]). That this indeed corresponds 
to Clement & Rowe s OoII component can be seen by the fact 


that this line is significantly displaced towards longer periods, 
compared to the Ool reference line (also shown in Figure[2]». 

As can be seen from Figure [3] there is a fairly tight relation 
between J- and Ks-band RRL amplitudes, which is described by 
the following expression: 


J[j = 2.6 ^ 2 . ( 2 ) 

Equation [2] was used to obtain the Ks-band equivalent of Equa¬ 
tion |T] OoII reference loci obtained in this way are overplotted 
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Fig. 2. Bailey diagrams for all the RRL stars known in the to Cen field, 
based on 7- {top) and A" s -band bottom light curves. The overplotted lines 
define representative loci of “OoII-like” {dashed lines ) and “Ool-like” 
{dotted lines) behavior for the ab-type RRL, obtained as described in the 
text. RRab stars with documented Blazhko effect are marked with open 
circles. Clearly, VI80 has too large an amplitude for its claimed RRc 
status; instead, we believe that the star is a contact binary of the W UMa 
type with twice the alternative “RRc-like” period. V151 is clearly an 
RRc star. 


on both panels of Figure [2] as dashed lines, clearly showing the 
expected period shift with respect to the (minority) “Ool-like” 
RRL stars in to Cen. 

Unfortunately, there are few stars in to Cen that belong to 
the Ool component (Clement & Row e[2000] >, and thus we can¬ 
not define the Ool reference line reliably based solely on our 
data. Therefore, to obtain an Ool reference locus in the NIR, 
we have used as a reference the V-band Ool line from ICacciaril 
|et al.| ( |2005| ) and |Zorotovic et al.| ( |2010] >, transformed into the 
NIR bandpasses by means of the following relations; 

.VL = 0.46,Vty, (3) 


Jl Ks =0.32yt“ /3 , (4) 

both of which are based on to Cen RRab stars, with visual am¬ 
plitudes coming from Kaluzny et al.](|2004[). The corresponding 
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Fig. 3. Relation between Ks- and 7-band amplitudes {top) and Ks- and 
V-band amplitudes ( bottom ) for to Cen RRab stars. The dashed lines are 
the fit to the data (eq.[2]an d0- 


J- and /ifs-band Ool lines are shown in Figure [2] as dotted lines. 
The presence of Ool-like RRL stars in to Cen is confirmed, and 
is especially evident in the 7-band Bailey diagram. 

These interrelations involving W, VI/, and VI Ks , compared 


to those previously presented by Feast et al. (2008), provide a de¬ 
scription of the amplitudes that is physically better motivated, as 
they more properly describe their asymptotic behavior towards 
vanishing light-curve amplitudes. 

It is still a matter of debate whether the position of stars on 
the Bailey diagram, and therefore the Oosterhoff type, is de¬ 
termined primarily by the metallicity of the cluster or by its 
HB morphology ( Contreras et al.|2005 2010 Clement & Rowe 


|2000| |Rey et al.||2000| and references therein). In this respect. 
u> Cen provides an interesting test of the relative importance of 
HB morphology and [Fe/H] in defining the Oosterhoff classifica¬ 
tion, given that the cluster possesses a predominantly blue HB, 
as in the case of most OoII clusters, but it also shows at least 
two well-defined metallicity peaks, separated by ~ 0.5 dex. To 
further investigate the possible impact of metallicity upon the 
distribution of RRL stars in the NIR Bailey diagrams. Figure [4] 
again shows the Bailey diagrams in both 7 and Ks, but separating 
the RRL into three metallicity bins, based on [Fe/H] values from 
Sollima et al. ( 2006] ) and Rey et al. (2000 1 . The bin sizes were 


chosen so as to ensure an equal number of RRab stars in each 
bin. For stars with [Fe/H] measurements in both catalogues, the 
metallicities of |Sollima et al.| were preferred. 

As shown in Figure [4] to Cen can still (not unexpectedly) 
be safely classified as predominantly an OoII GC, which is in 
perfect agreement with its predominantly blue HB morphology. 
This figure also shows that the NIR Bailey diagrams for to Cen 
ab-type RRL stars is largely insensitive to metallicity, in the 
sense that the average OoII and Ool loci, defined as above, ad¬ 
equately represent the mean loci occupied by RRL in each in¬ 
dividual metallicity bin also. It is interesting to note, however, 
that the Ool/OoII number ratio does change with metallicity: 
in the lowest-metallicity bin, there are virtually no RRL stars 
close to the Ool locus; conversely, in the highest-metallicity bin, 
RRab stars are more evenly divided between OoII and Ool sub¬ 


classes. This largely confirms previous results from Rey et al. 


(|2000]> based on visual data, even though the impact of metallic- 
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Fig. 4. As in Figure[2] but separating the RRL stars into three metallicity bins, from lowest (top) to highest (bottom). The sizes of the bins were 
chosen so as to include the same number of RRab stars in each. 


ity at low [Fe/H] was stronger in their case. Thus, HB morphol¬ 
ogy and metallicity appear both to be important in determining 
Oosterhoff type. 

Interestingly, Figure 9 in Rey et al. |2000] l suggests that 
the to Cen ab-type RR Lyrae with the longest periods (i.e., 
P > 0.71 d) all belong to the metal-intermediate bin, but this is 
not immediately obvious from our plot. We tested whether this 
difference is simply due to the different metallicity bins selected 
by Rey et alT] and our study. Indeed, according to our analysis. 


even using the Sollima et al. (20061 metallicities, the long-period 
RRab’s in to Cen RRab’s do indeed occupy a very limited range 


in [Fe/H], between -1.81 and -1.64. This can be understood in 
terms of the evolutionary scenario proposed by |Lee| ( 1991) , ac¬ 
cording to which the HB component over this metallicity range 
is extremely blue, thus preferentially producing overluminous, 
long-period RRL stars. 

3.2. Low-amplitude RRc stars 

Some RRc stars in to Cen are known to have very low amplitudes 
in the visible, at the level of fHy <0.1 mag ( |Kaluzny et al. 2004) . 
Specific examples include V281 (3\y ~ 0.07 mag), V344 (Jtv at 
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Fig. 5. Finding charts for the 4 new RRL candidates that were discovered in our VISTA images. North is up and East to the left, and each stamp 
covers a field of 1 arcmin 2 . The discovered candidates are the stars at the center of each stamp. The first two stamps correspond to 2 new RRL 
members of to Cen, based on their position in the CMD (Fig. [TJ, whereas the remaining ones, according to the same CMD, are field interlopers, in 
the cluster background. 

Table 2. New RR Lyrae candidate stars'"’ 


ID 

RA 

(J2000) 

DEC 

(J2000) 

(J) 

(mag) 

(K s ) 

(mag) 

P 

(days) 

Type 

d 

(arcmin) 

Memb. 

status 

NV455 

13:27:53.94 

-46:55:43.93 

13.132 

12.760 

0.932469 

RRab 

34.92 

m 

NV456 

13:22:14.49 

-47:24:21.64 

13.530 

13.270 

0.383516 

RRc 

46.33 

m 

NV457 

13:29:54.56 

-47:50:46.00 

15.530 

15.265 

0.508593 

RRab 

38.45 

f 

NV458 

13:30:00.09 

-47:13:05.63 

14.772 

14.441 

0.620313 

RRab 

36.23 

f 


Notes. {o) Columns have the same meaning as in Table[I] 


0.08 mag), V357 (Jt v as 0.05 mag), and V399 (LTlv ~ 0.05 mag). 
We do not expect to recover the variability of these stars from our 
NIR observations, because the amplitudes in the NIR are even 
lower than in the optical (Smith 1995[ >. Indeed, none of these 
stars were found to be variable in our data. 

In this context, we also note that clean light curves for V339 
and V340 could not be recovered either, even though their V 
amplitudes (Jt v as (). 11 mag and 0.17 mag, respectively) are per¬ 
haps high enough that one might have expected non-negligible 
J and K$ amplitudes. A possible explanation for their non de¬ 
tection in our data is their marked multiperiodic behavior, which 


renders even their V -hand light curves somewhat noisy (Kaluzny 
|et al.|2004) . 


3.3. New RR Lyrae stars 

To detect additional, previously unknown variable stars in and 
around to Cen, we performed a new search for variable stars over 
the complete area covered by our observations. Our procedure 
consisted in considering as candidate variables all those stars 
whose light curves implied root-mean-square magnitude devi¬ 
ations that were higher than the average trend for non-variable 
stars at the same magnitude. The search was independently per¬ 
formed in J and K$. With the goal of specifically detecting RRL 
and other periodic variables, the four most prominent periods for 
each star, as returned by the ANOVA algorithm (Schw arzenberg-| 
|Czemy|[1989| l, were used to plot phase-folded variables. From 
this search, four variable star candidates have periods and ampli¬ 
tudes in accordance with their being RRL stars. Figure [5] shows 
the finding chart for these new discovered RRL candidates. None 
of them appears to present the Blazhko effect. The remaining 
candidates, which may belong to different variability classes, 
will be discussed in future papers of this series. 

Table |2]gives the ID (following the numbering scheme from 
the C13 online catalogue), coordinates, intensity-weighted mean 
J and K s magnitudes, periods, type, and the distance to the cen¬ 
ter of the cluster. According to their position in the CMD, NV455 


and NV456 are likely members of the cluster, whereas the other 
two stars deviate largely from the main HB locus - and since 
they are fainter, they are almost certainly field stars in the cluster 
background. J- and /Ts-band light curves are presented in Fig¬ 
ure [6] According to their light curve shapes, in addition to their 
periods (Table [2} and positions in the Bailey diagram (Fig. [2}, 
we assign an RRc type to the new variable NV456, and an RRab 
type to variables NV455, NV457, and NV458. As it happens, 
NV455, with a period of about 0.932 d, is the RRab star with the 
third longest period known to date in to Cen (see Fig.[2ji. 

Recently, Fernandez Trincado et al. < 2013] ) also performed a 
variability search in the V -band over an area of 50 deg 2 around 
to Cen, with the main aim to detect RRL stars and in this way find 
evidence for tidal debris from the putative dwarf galaxy progen¬ 
itor of the cluster (e.g., Dinescu|2002||Altmann et al.|2005 |Car-| 


retta et al.|2010 Valcarce & Catel an|201 1[ >. They found 29 new 

RRL stars, and recovered 7 previously known ones. Of these, a 
total of 11 are in the area covered by our observations, and they 
were all independently recovered in our study as well. Their re¬ 
maining RRL are all outside to Cen’s tidal radius, and they are 
reported further on in |Fernandez-Trincado et al.| ( )20T4| . 


4. Summary 


Based on previous to Cen variability studies, from Bailey (1902) 
until |Kaluzny et ah] ( |2004| >, C13 report that there are 193 RRL 
stars in the field of the cluster. Among these, one finds 92 RRab, 
100 RRc, and 1 star classified as RR? (V80). None of the 5 new 
RRL claimed by [Weldrake et al.[ ( j2007| ) were included in the 
Cl3, catalogue, because, as we have seen, only two of them are 
really new discovered - and for these two stars, ID-91 and ID-99, 
an RRL classification still remains uncertain. 

On the other hand, ID-74, which [Weldrake et al.| ( ]2007] ) clas¬ 
sified simply as “pulsating,” was included in the C13 catalogue 
as a possible RRL, based on its period (0.6671 d). In this way, it 
was also given the official ID V433 by C13. However, we have 
been unable to recover the star’s variability using our NIR data. 
Moreover, its NIR magnitude and color place the star close to but 
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Fig. 6. J and A's-band light curves of the newly discovered RRL stars that are likely cluster members. 


to the red of the base of the red giant branch, rather than on the 
classical IS. Recently, Navarrete et al. ( 2013|) found a misiden- 
tified RRab star, V411, in the Weldrake et al. catalogue. There¬ 
fore, prior to this publication, 194 RRL stars in all (93 and 100 
RRab and RRc, respectively) have been listed for the co Cen field, 
where we have not included any of the “new” candidates in the 
|Weldrake et al.[ catalogue in this number count. 

Considering the information obtained with the VISTA obser¬ 
vations, the light curves of some RRLs were recovered for the 
first time, leading to changes in the variability classification. In 
particular, we have been able to provide the first light curves ever 
for V80, thus establishing that it is an RRc star. V151, V173, and 
V177, all of which had previously been considered ab-types, are 
shown to be RRc stars, according to our data. On the other hand, 
we have found that V180 is not an RRc, as previously thought, 
but rather a contact binary of the W UMa type. Furthermore, 
4 new variables were found using VISTA: 3 RRab and 1 RRc 
stars. Therefore, there are now 197 known RRLs in the co Cen 
field, including 93 RRab and 104 RRc stars. 

From these 197 RRLs, there are 4 RRab (V181, V283, 
NV457, and NV458) and 2 RRc (V168 and V183) stars that have 
magnitudes and colors consistent with being field interlopers. In 
addition, if we exclude V68 and V84 because of their controver¬ 
sial classification (with the distinct possibility that they may be 
ACEP stars, rather than RRLs), the final census of likely RRL 
members of co Cen includes 88 RRab and 101 RRc stars, for a 
grand total of 189. 

It should be noted that it was not possible for our team to ob¬ 
serve some distant RRLs such as V171, V178, V179, and V182, 
because they fall outside the VIRCAM@VISTA FoV; however, 
they are located outside the tidal radius of the cluster, and are 
therefore most probably field stars. It is expected that these num¬ 
bers will increase further, once the new variables found by Fer- 
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jnandez Trincado et al.| ( j2013| ) in the co Cen neighborhood are 
published. 

In a forthcoming paper, we will present a detailed analy¬ 
sis of the NIR period-luminosity relations followed by RRLs as 
well as type II Cepheids and SX Phoenicis stars in the co Cen 
field, based on our VISTA observations. This will be followed 
by a paper providing the complete, updated catalogue of variable 
stars in co Cen, including the full list of RRL, type II Cepheids, 
SX Phoenicis stars, eclipsing binaries, and anomalous Cepheids, 
along with their NIR properties. 
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Appendix A: Notes on individual RR Lyrae stars 


V19 and V21: In the C13 catalogue, these variables are listed 
as possible non-members, based on their anomalous Fourier pa¬ 
rameters compared with other variables with similar periods 
( [Clement & Rowe|2000] l. However, and as indicated in Table [1] 
the proper motion studies of vLOO and B09 found a membership 
probability P M of (100%, 100%) and (100%, 99%) for V19 and 
V21, respectively. The membership status for these two stars is 
further supported by their position in our NIR CMD, which lo¬ 
cates them in the middle of the RRL IS (see the bottom panel of 

Fig-0- 

V52: This star appears to be much brighter than other RRab 
stars in u> Cen, as can be noted in Figure 0 Its light curve has 
high photometric errors despite its high luminosity, thus suggest¬ 
ing contamination due to blending. Indeed, vLOO reported that 
V52 may have an unresolved companion, which would account 
for the fact that it is brighter and has a lower amplitude than other 
RRL variables with the same period. According to vLOO and 
B09, the membership probability of this star is 42% and 53%, 
respectively. Based on its position in our NIR CMD, the pro¬ 
jected distance to the cluster center, d ~ 2', and the high amount 
of noise in the light curve, V52 is most likely a member of the 
cluster in a very marked blend. Inspection of Hubble Space Tele¬ 
scope (HST) images ( [Anderson & van der Marel|2010| confirms 
the presence of a companion that is unresolved in our images, 
located just about 0.5"from the RRL star, but brighter by about 
0.5 mag in V, and bluer by about 0.25 mag in B - V. 

V56 and V168: The field status of these two stars was based 
on a radial velocity study by Liller & Tokarz ( 1981j ). Subse¬ 
quently, V56 was again considered a member of the cluster based 
on its high membership probability of 98% and 100% found 
by vLOO and B09, respectively - a status in excellent agree¬ 
ment with its position in our NIR CMD. These inconsistent con¬ 
clusions can probably be explained as follows. When assessing 
membership status of an RRL variable from its radial velocity, it 
is necessary to take into account the phase in its light and veloc¬ 
ity cycle. [Liller & Tokarz| obtained this information from R. I. 
Dickens (1980, priv. comm.), based on an epoch established in 
1974. The observations of lLiller & Tokarz! were made in 1978. 
A period change study by [Jurcsik et al. | ( j2001 [ > indicates that V56 
undergoes random period changes. Thus it is highly likely that 
the light curve phase assumed by |Liller & Tokarz] was incorrect. 
On the contrary, all proper motion studies agree in classifying 
VI68 as a non cluster member, which also agrees with our pho¬ 
tometric classification as a field star: Figure 0 clearly shows that 
this RRL star is a background interloper. 

V68 and V84: These two stars have controversial mode 


classification in the literature. First thought as RRLs by Bailey 
( |1902|l, they wer e later suggested to belong to the ACEP group 
by Nemec et al.j( |1994[ ), because of its brighter optical magni¬ 
tudes compared to the other RRLs of the cluster. ACEP stars 
are Population II stars with periods between 0.4 and 2.4 days, 
that do not follow either the RRL or the type II Cepheid PL 
relations. Most of them are found in dwarf spheroidal galaxies 
( Nemec et al.|1 994) and the Large Magellanic Cloud (Soszynski 
|et al,|2008[>, a few ACEPs are known in Galactic globular clus 


ters (s ee |Zinn & Dahn|197'6j Sollima et al,|2010||Qsbom et al.| 


2012 |Arellano Ferro et al.||201 3|> an d only one in the Galactic 

field (XZ Ceti; |Teays & Simon|1985j|Szabados et al.|2007| l. 

V68 appears as an RRc cluster member star according to its 
proper motions and its position in the CMD (see Table 0 and 
the bottom panel of Fig. 0. Although it has the longest period 
among the RRc stars of the cluster (0.56 days), which is longer 
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than even the short-period end of the RRab stars’ period distri¬ 
bution, the amplitude of its light curve is small enough to place 
it within the RRc stars group (see Fig. [2]i. V68 also follows the 
PL relations for RRc stars in J and Ks (Navarrete et al. 2014b, 
in preparation). This star, on the other hand, also appears to fol¬ 
low the optical (B, V) PL relations for ACEP stars, as derived by 
|Nemec et al.| ( |1994] >. 

On the other hand, V84 was discarded as a cluster member 
by vLOO, with a membership probability = 0%. Unfortu¬ 
nately, we were unable to confirm this result using proper mo¬ 
tions from B09, since V84 is outside their FoV. Our NIR CMD 
shows that V84 is well located in the RRL instability strip, open¬ 
ing the possibility that it is a cluster member. However, as |Del| 
|Principe et al.[ ( j2U06 ) pointed out, V84 does not follow the NIR 
PL relation for RRab stars. This result is confirmed using our 
PL relations (Navarrete et al. 2014b, in preparation), with V84 
appearing brighter than other stars with the same period. Con¬ 
versely, Nemec et al.]( 19941 found that V84 was placed slightly 
below the B- and V -hand PL relations for ACEPs, which could 
be explained if the star was not an ACEP but instead an RRL star 
leaving the HB. 

Using our NIR observations, we cannot adopt a definitive 
classification for either of these stars. If they are RRLs, the most 
probable scenario is that V68 is an RRc member and V84 is an 
RRab star from the field. However, the ACEP classification can¬ 
not be completely discarded. Indeed, V68 and V84 follow the 
/(Vband PL relation for ACEPs in the Large Magellanic Cloud 
derived by jRipepi et al. ( |2014) >, assuming a distance modulus for 
to Cen of (m - M)o = 13.63 + 0.1 mag - which is in good agree¬ 
ment with the one derived by Matsunaga et al. ( 2006j l using type 
II Cepheids. If these two stars are indeed ACEPs, they should be 
members of the cluster. 

V80: This star was detected in the first variability study on 
to Cen carried out by Bailey ( [1902 1. In that work, V80 was re¬ 
ported as difficult to measure because it is distant from compari¬ 
son stars, and accordingly only sparse measurements of differen¬ 
tial magnitudes were published, without an accompanying light 
Bailey| suggested a period of 0.45 or 0.31 days for this 


curve, 
star, w 


terewith it was tentatively classified as an RRL-type star. 
Based on our wide-field VISTA observations, we determined its 
period, P = 0.37718 d (Table[I]i, and derived its complete, phase- 
folded light curve (Fig.|A.l|. 


more than 10 stars inside a 1" radius around its coordinates. A 
similar situation occurs with VI39, having a neighbor located at 
only 0.432", according to the Anderson & van der Marel cata¬ 
logue - which is a difference of less than 2 pixels, given the VIR- 
CAM pixel scale. Hence, this star is almost certainly brighter in 
our photometry than it would be, in case it were isolated. 

V135 is suffering the influence of a disturbing neighbor, lo¬ 
cated 2.5" away, which is saturated. This situation was previ¬ 
ously reported by Samus et al. ( |2009| ), with both stars being la¬ 
beled in that study as a visual binary. 

V143: It is certainly a member of the cluster according 
to vLOO, with 99% membership probability. B09 give instead 
Pf, = 90% for this star. However, it appears a few tenths of a 
magnitude brighter than other RRab stars with the same color. 
A visual inspection of this star in some of our images suggested 
to us that its anomalous position in the CMD is not due to a 
blend effect of an unresolved system, but rather to a saturated 
companion located close by (i.e., at a distance of ~ 3.4") that 
is somehow disturbing the photometric quality in its neighbor¬ 
hood. Indeed, this situation was explicitly reported by |Shokin| 


According to its coordinates, ID-53, P = 0.377 d, from Wel- 
[drake et aT7[ ( j2007) > appears as the nearest star to V80 among the 
list of new variables discovered by those authors. However,|Wel- 
drake et al. did not associated ID-53 to V80, probably because a 


period value had never previously been published for V80 - and 
accordingly, they claimed that their ID-53 was a newly discov¬ 
ered variable star. Instead, our independently derived period and 
coordinates alike indicate that V80 and ID-53 are one and the 
same star. 

V80 is not present in any proper motion study. Its derived 
period and position in the CMD indicates that V80 is an RRc 
star that is also a member of the cluster. 

VI18, V135, and V139: These three stars appear well lo¬ 
cated within the RRab region in the IS (see Fig.jTj). Based on the 
vLOO and B09 proper motion studies, VI18 has a membership 
probability of P^ = (100%, 96%), respectively, whereas VI39 
has P M = (100%, 90%). VI35 is only present in B09, with a 
membership probability of 97%. However, all appear brighter 
than the other RRab stars with similar periods in the J and Ks 
PL relations (Navarrete et al. 2014b, in prep.). 

VI18 is located in a very crowded region. According to the 
HST catalogue of Anderson & van der Marel (|2010[>, there are 


& Samus’ ( jl996] >, where the presence of a bright star located 
southwest from the position of V143 is remarked upon. There¬ 
fore, despite its anomalous position in the CMD, we suggest that 
V143 is indeed a cluster member. 

V151: Listed as an RRab star by Kaluzny et al. ( [2004 1 with 
a period of 0.4078 days, this star would appear to be the shortest- 
period fundamental-mode RRL in to Cen. However, despite the 
fact that our estimated period (P = 0.407756 d) is in excellent 
agreement with the one derived by Martin ( 1938) , we notice that 
its light curve shows a very small amplitude, which is among 
the lowest for RRab stars in our study. Small amplitudes are un¬ 
usual for short-period RRab stars. Figure [2] shows the position 
of V151 in the Bailey (period-amplitude) diagram, as compared 
to other to Cen RRL stars. This plot clearly shows that V151 is 
well located among the RRc-type variables. In addition, its posi¬ 
tion in the PL relation (Navarrete et al. 2014b, in preparation) is 
also consistent with first-overtone pulsation. Thus, we conclude 
that this star should be classified as an RRc. We note that lMartml 
( |1938[ l also classified V151 as an RRc star, hence the erroneous 
RRab type adopted by |Kaluzny et al. is likely a typo, which was 
later inadvertently reproduced in the Cl3 online catalogue. V151 
was not included in the proper motion studies of vLOO and B09; 
however, our NIR CMD suggests that this variable likely belongs 
to the cluster. 

VI72, V173, V175, V177, and V180: |Wilkens[ ( jT965] > an¬ 
nounced the discovery of V171-V180. Of these 10 stars, V176 
was subsequently discarded as a non-variable star, V174 clas¬ 
sified as an EB, and the remaining stars classified as RRab’s, 
except for VI80, which was classified as an RRc instead. To 
the best of our knowledge, no further variability studies of these 
stars have been carried out, except for V172 and V177 which 


were included in the observations by Weldrake et al. (20071 as 


ID-18 and ID-49, respectively. None of them has had a period 
assigned in the C13 catalogue. 

Four of Wilkens’s stars, namely, V171, V175 (without a pre¬ 
vious type assigned), V178, and V179 are all very far from the 
cluster center, and thus were not observed by us, despite VIR- 
CAM’s large FoV. Consequently, only light curves for V172, 
V173, V177, and V180 were derived, and they are shown in Fig- 
Only VI80 has membership probability derived based 


A. 1 


on proper motions studies, being considered a non-member of 
the cluster with a null membership probability by vLOO. Accord¬ 
ing to their positions in the CMD (bottom panel of Fig. |T]», we 
suggest that V172, V173, and V177 are cluster members. We 
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Fig. A.l. VISTA NIR light curves for all the RRLs that lacked a derived light curve or period in the literature. For each star, the upper panel shows 
the 7-band curve, whereas the lower one displays the K$ band instead. Our derived period for each star is provided in each panel. Note that the 
V180 light curve is consistent with a W UMa classification. Variables ID-91 and (especially) ID-99, from |Weldrake et aT7] ( |2007[ >, have light curve 
shapes that are different from all the other RRab in a> Cen, and are suspected field stars. 
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confirm the ab-type classification for V172, but at odds with 
Wilkens, we did not find an RRab behavior for either V173 or 
V177, and we suggest instead an RRc classification for both 
these stars. This is consistent with their positions in the Bailey 
diagram, as can be seen from Figure |2| V175, according to A. 
K. Vivas (2014, priv. comm.; see also Fernandez Trincado et al. 
|2013j ) is a foreground c-type RRL. As can be seen in Figure [T] 
VI80 is located in the RRab region, redder than any other RRc 
star. Moreover, it has a large amplitude compared to other RRc 
stars (Fig.|2|. Our estimated period produces a light curve clearly 
showing 2 different minima. Accordingly, we propose to classify 
VI80 as a contact EB of the W Ursae Majoris type, rather than 
as an RRc, which is consistent with its fairly red color and the 
fact that this star is not a cluster member. 

V181 and VI83: V181, VI82, and VI83 were discovered 
by Wesselink (1969, unpublished, priv. comm, to H. Sawyer- 
Hogg), who indicated the stars’ coordinates and periods, but 
without providing light curves or variability types for these stars. 
Clement et al. ( 2001) >, based on Wesselink’s periods, classified 


Therefore, V349 and V351 could not be recovered using PSF 
photometry in our VISTA images. 


Appendix B: RR Lyrae stars from Weldrake et al. 
12007) 


V181 and VI82 as RRab’s, and VI83 as an RRc. These stars 
were outside the FoV in subsequent investigations by other au¬ 
thors. Unfortunately, VI82, at 59.7 arcmin from the center of the 
cluster, is also outside the FoV of VISTA. 


Figure A.l shows the first complete light curves for V181 Both stars’ light curves are shown in Figure A.l 


Weldrake et al.. f2007| reported the discovery of five new RRL 
stars in to Cen: ID-53, ID-91, ID-99, ID-144, and ID-145. By 
carefully checking these stars’ positions and periods, we have 
been able to establish that three out of these five variables were 
not really new discoveries. More specifically, according to our 
analysis, we were able to establish that ID-53 = V80, ID-144 = 
V272, and ID-145 = V15. 

ID-91 and ID-99, on the other hand, do appear to be bona- 
fide new variables. Weldrake et al. (2007) listed optical mag¬ 
nitudes of V ~ 17.1 and 17.45 mag for these two stars, re¬ 
spectively, which is consistent with their being field interlopers. 
Unfortunately, neither star appears in the proper motion cata¬ 
logues of vLOO and B09. As far as their variability status is con¬ 
cerned, Weldrake et al. suggested that ID-99 could be either a 
long-period RRc with a 0.671-d period or an EB with a period 
twice as long. For ID-91, in turn, they assigned an RRab type. 

according to 


and V183. They are ab- and c-type RRL stars, respectively, 
in agreement with the periods originally derived by Wesselink. 
These two stars are very distant from the cluster center (though 
still inside its tidal radius), and thus they are absent from proper 
motion studies of the cluster. However, based on their positions 
in the CMD (see the top panel of Fig. [TJ, they both appear to be 
field RRL stars. 

V268: vLOO claimed that V268 is a field star with a member¬ 
ship probability of 0%. However, we found some inconsistencies 
in this membership status using our coordinates and those from 
the literature: in fact, the nearest star to V268 in vLOO’s cata¬ 
logue has a membership probability of 70%. Moreover accord¬ 
ing to B09, this star has a 99% membership probability. Using 
the NIR CMD to figure out the membership status of this star, 
it is found that V268 is located well inside the RRL instability 
strip, as the other RRab members of the cluster (Fig. |TJ bottom 
panel). We thus conclude that V268 is indeed a cluster member. 

V283: This star was studied by neither vLOO or B09. Based 
on its faint visible magnitude, Kaluzny et al. ( |1997 1 claimed that 
it is a background halo RRL. This field status is confirmed by 
our photometry, which locates the star a few magnitudes below 
the HB, even fainter than the turn-off point level (see CMD in 
the top panel of Fig.|T). 

V349 and V351: These two RRc variables were discovered 
by Kaluzny et al.] ( j2004] l and are located ~ 1.5 arcsec away 
from the cluster center. The huge amount of crowding at their 
positions did not allow those authors to present properly cali¬ 
brated light curves. Interestingly, Kaluzny et al. noted that V351 
is a multiperiodic variable, and suggested that it may show non¬ 


which the RRab status of ID-91 seems reasonable, whereas the 
RRL nature of ID-99 remains a possibility. However, if indeed 
an RRL, ID-99’s long period might more naturally suggest an 
RRab classification, as opposed to the RRc type suggested by 
|Weldrake et al.] To further check their classification, ID-91 and 
ID-99 were plotted in the Bailey diagram (see Fig.|2|. Their lo¬ 
cation suggest that both stars are either extreme examples of Ool 
RRab stars, and/or show a modulation in their light curves (e.g., 
the Blazhko effect), and/or are unresolved blends. The latter pos¬ 
sibility receives some support from the positions of the stars in 
the CMD (Fig. |T]), which show that both stars are slightly redder 
than expected for RRL stars. However, these stars are also quite 
far away (~ 25') from the cluster center, and inspection of our 
images reveals that they are in relatively unpopulated areas, thus 
making the possibility of blends seem unlikely. Unfortunately, 
HST images of these fields are not available, and thus we cannot 
conclusively establish that unresolved companions to these stars 
are indeed absent. In like vein, the presence of the Blazhko ef¬ 
fect is not supported by either our photometry (Fig. |A. 1 1 ) or the 
|Weldrake et al.| ( |2007| photometry. 

Based on the astrometry that is provided along with the 
VISTA images (which is tied to the 2MASS system), we found 
some erroneous matches between some of the known variables 
in the catalogues of Weldrake et al.|(|2007) and Kaluzny et al. 
( |2004l i: ID-115 = V2 64 (instead of V 48) and ID-135 = V266 
(not V356). Moreover, Weldrake et al. erroneously claimed that 
ID-133 was the same star as V144 from the IClement et al.l 
P001 catalogue. However, as already reported by Navarrete 


radial pulsation modes, as previously proposed by Clement & 
|Rowe POOP) for several other RRL stars in to Cen. Using the 


pf al.|(2013) , ID -133 is indeed a different variable star. C13 has 


coordinates and periods provided by |Kaluzny et al.| ( |200 4), these 
two variables were not found neither in our ALLFRAME nor in 
our DoPhot photometry. The matching radius used was 1 arcsec. 
Increasing the matching radius until 10 arcsec still did not return 
any variable stars. The last check step was to consider all the 
stars within 30 arcsec around the listed coordinates, and check 
for the presence of periodicity in this enlarged sample using the 


Analysis of Variance statistic (ANOVA; Schwarzenberg-Czerny 


11989). Again, no periodic variables were detected in this way. 


already added this new RRab star in her online catalogue, where 
it is now officially listed as V411. 

C13 has included ID-74 of Weldrake et al.|p007), a possibl e 
RRL based on its period (0.6671 d, according to |Weldrake et al.) , 
as V433 in her catalogue. This star has V ~ 16.5 and (V - I) » 
1.0 mag, which suggests that it is a field star. Indeed, according 
to our NIR magnitudes, V433 star is placed above the cluster’s 
turn-off point, above the main sequence and close to (but redder 
than) the base of the red giant branch (see the top panel of FigJTj). 
Still, the star does appear to show some short-period variability, 
with a period of (0.6681 + 0.0004) d, and an amplitude that is 
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likely not higher than 0.05 mag in Ks, again preventing us from 
classifying it as an RRL. The star’s red color is unlikely to be 
caused by a blend with a red star, since it is ~ 30' away from 
the cluster center, and without any obvious close companions 
in our images. According to Da Costa & Coleman (2008), the 
star is a radial velocity non-member. In conclusion, we follow 
|Weldrake et al.| ( [2007[ ), and consider V433 as a variable star with 
an uncertain variability type. 
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